Area coverage of radial Levy flights with periodic boundary conditions 
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We consider the time evolution of two-dimensional Levy flights in a finite area with periodic 
boundary conditions. From simulations we show that the fractal path dimension df and thus the 
degree of area coverage grows in time until it reaches the saturation value d/ = 2 at sufficiently 
long times. We also investigate the time evolution of the probability density function and associated 
moments in these boundary conditions. Finally we consider the mean first passage time as function 
of the stable index. Our findings are of interest to assess the ergodic behavior of Levy flights, to 
estimate their efficiency as stochastic search mechanisms and to discriminate them from other types 
of search processes. 

PACS numbers: 05.40.-a,02.50.-r 



I. INTRODUCTION 



Levy flights are random walk processes, in which the 
lengths of individual jumps are distributed according to 
a probability density \{x) of the asymptotic power-law 
form 0-3 



\{x) 



\l+a ■ 



< a < 2, 



(1) 



where er is a scaling factor of physical dimension length, 
[<x] = cm. The resulting motion is spatially scale- 
free due to the divergence of the jump length variance, 
/ x 2 \(x)dx. Levy flights were popularized by Benoit 
Mandelbrot, who named them after his teacher, French 
mathematician Paul Levy (BJ . They have been mainly ap- 
plied in the modeling of search processes, following the 
original idea by Klafter and Shlcsinger [fj: while regu- 
lar random walks in one or two dimensions have a high 
probability to return to already visited sites, Levy flights 
combine thorough local search with occasional long ex- 
cursions, leading them to areas which likely have not 
been visited before. This strategy reduces unnecessary 
oversampling and thus represents an advantage to the 
searcher. Indeed, trajectories consistent with the long- 
tailed jump length pattern ([1]) have been reported for 
various animal species [7H1G] . Even for molecular search, 
Levy flights may emerge from the topology of the search 
space UJ. Due to a mix of modern means of trans- 
portation also human motion behavior is characterized 
by travel lengths with power-law distribution [l2| ■ Long- 
tailed distributions of relocation lengths also change sig- 
nificantly the distribution pattern of diseases, as regular 
diffusion fronts arc broken by, for instance, long-distance 
air travel, thus carrying the disease to completely decou- 
pled places [HI]. 

For land-based animals typical search processes are es- 
sentially two-dimensional. In many cases we may also 
neglect the vertical dimension for many birds or fishes, 
when the lateral extension is considerably larger than 



the maximal height/depth difference of the trajectory. 
In such a two-dimensional search space, given its fractal 
dimension df = a a single Levy flight trajectory can- 
not fully cover the search space in absence of boundaries. 
Thus, some small area element in the search space may 
be hit by one single trajectory while it may be missed 
by another. In many cases, however, the search space is 
bounded. For instance, animals only search for food in 
their own territory, or the business travel patterns of an 
individual are confined to a certain country or continent. 
Moreover it is often relevant to have information about 
the area coverage of a Levy flight process at finite times. 
For instance, a predator often does not cover its entire 
territory on a single day. It is therefore of interest to 
explore the time evolution of the area coverage of Levy 
flights: how long does it take for the animal to efficiently 
explore its entire patch, or a disease to reach every little 
town in a country? How large is the area coverage at a 
given time? These questions of area coverage are directly 
connected with the ergodic properties of Levy flights. 

Incomplete area coverage of Levy flights is also of in- 
terest in view of their ergodic properties. Ergodicity and 
its violation have recently received considerable inter- 
est, following the possibility to record single trajectories 
of molecules or small tracers in biological matter [lij ]. 
Thus, for anomalous diffusion processes with underlying 
scale-free motion a weak ergodicity breaking is observed 
fl~5l HH, such that the mean squared displacement ob- 
tained from the time average over a single trajectory is a 
random variable and significantly different from the cor- 
responding ensemble average, as shown for the subdiffu- 
sive continuous time random walk process [l~7l - [2l| . This 
behavior contrasts the ergodic behavior of other anoma- 
lous diffusion processes such as fractional Brownian mo- 
tion [U [H[ . While weak ergodicity violation occurs for 
processes with diverging characteristic time scales, ergod- 
icity is also violated for Levy flights in infinite media, due 
to the incomplete area coverage even for infinitely long 
trajectories [24[. Confined Levy flights, in contrast, are 
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FIG. 1: (Color online) Trajectories of radial two-dimensional 
Levy flights for different values of the stable exponent a, 
showing 400 steps each in a 1000 x 1000 box under peri- 
odic boundary conditions. Starting point was the center of 
the square area. Note the significantly longer excursions for 
smaller values of a. We chose a = v2. 



crgodic [25 1 . 

In this paper, we consider radial Levy flights in two- 
dimensional square areas with periodic boundary condi- 
tions. We investigate the time evolution of the effective 
fractal dimension of the flights for varying stable index 
a. As a complementary measure for the area coverage we 
analyze the time evolution of the mean squared displace- 
ment up to its saturation. Moreover we analyze the time 
evolution of the probability density function and its mo- 
ments. We also pursue the question on the typical time 
for the Levy flight to first reach or pass the boundary of 
the square interval. 

The paper is organized as follows. In section II, 
we briefly review the mathematical foundation of Levy 
flights. In section III, we consider radial two-dimensional 
Levy flights in square boxes and numerically study their 
area coverage in terms of their fractal dimension, their 
probability density function, evolution of moments, and 
mean first passage times. Section IV summarizes the re- 
sults. 



II. LEVY FLIGHTS 



For jump length distributions of the type ([TJ, the vari- 
ance (x 2 ) = f_ x 2 \{x)dx diverges, while fractional mo- 
ments (\x\ s ) of order < 6 < a exist This property 
directly carries over to the random motion itself and is 
responsible for several peculiar phenomena. Due to this 
lack of a finite variance Levy flights have a fractal graph 
dimension df = a |2| . This property causes the sample 
path of a Levy flight to be characterized by occasional 



long jumps, connecting clusters of more localized jumps. 
Unvisited holes exist on all scales within the sample path. 
This scale- free behavior gives rise to the fact that clusters 
host smaller clusters, as shown in Fig. [T] Moreover, the 
discontinuous jumps on arbitrarily large scales induce a 
principal discrepancy between first passage and first ar- 
rival events. For instance, while their first passage behav- 
ior on a semi-infinite domain obeys the Sparre Anderson 
universality [26l . |27| , their first arrival behavior is signif- 
icantly reduced with decreasing a H|, as Levy flights 
strongly overshoot a point target [27 1 . On finite domains 
the first passage behavior is also modified [HI, [2t| . An- 
other interesting effect is the occurrence of multimodal 
distributions for Levy flights in the presence of steeper 
than harmonic potentials [3(| HH . 

Analytically, the statistics of the Levy flight are dic- 
tated by the detailed distribution of individual jump 
lengths. However, on large scales — when large search 
spaces are explored for sufficiently long periods of time — 
we obtain a universal description. Displacements on 
large time scales are characterized by Levy a-stable laws, 
which in the symmetric case are defined in terms of their 
characteristic function (exp(ifcx)) 0-[l, H3j 



\(k) 



\(x)e lk!X dx = exp(-a a \k\ a ). (2) 



Again, a is a typical length scale. From the details of 
the microscopic jump length distribution, only the tail 
property (fT]) is directly passed on to the asymptotic sta- 
ble law ^) of macroscopic displacements, including in 
particular the exact tail exponent a. Levy stable laws 
emerge as the limit of the sum of independent, identi- 
cally distributed random variables with diverging vari- 
ance, by virtue of the generalized central limit theorem 
0, H, H, HH . For < a < 2 this leads to the asymptotic 
power-law {J), while in the limit a = 2, we recover a 
Gaussian distribution for X(x) with finite moments of all 
orders. 

Here, we consider two-dimensional Levy flights in 
square boxes of edge length 2a, with periodic bound- 
ary conditions. This means that each time the particle 
crosses one edge of the box, it will enter from the oppo- 
site edge. All Levy flights start from the origin in the 
center of the square box. For each step in the two di- 
mensional random walk we draw a (signed, symmetric) 
flight distance r, with distribution A(r) as in @, and 
independently a flight direction 9, uniformly distributed 
within [0;7r]. The projection of displacements onto the 
x-axis is thus Ax = rcos(0). We call such processes 
radial two-dimensional Levy flights. Such radial Levy 
flights have been used to analyze search processes (see, 
e.g., Ref. @). We note that they are slightly different 
from two-dimensional Levy flights defined in terms of the 
two-dimensional characteristic function exp(— [cr*] Q |k| Q ) 

In what follows we calculate parameters such as the 
fractal dimension, the mean first passage time, the mean 
squared displacement, the probability density function 
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Input a Number of steps Output a 

0.7 10,000 0.68 ±0.06 

1.2 50,000 1.16 ±0.04 

1.7 10,000,000 |l.65±0X)3 

TABLE I: Value of the stability index a of free Levy flights. 
We list the input a used to generate the trajectory and the 
output value for a as determined from the box counting 
method. The number of steps is an indication for the number 
of steps of the simulated Levy flight necessary to determine 
a to reasonable accuracy. 

and its other moments and follow their temporal evolu- 
tion. 



III. RESULTS AND DISCUSSION 
A. Fractal dimension 

The fractal dimension df indicates how completely a 
fractal graph fills the available embedding space. For a 
Levy flight evolving in our finite area the value of the frac- 
tal dimension will initially be smaller than two. Starting 
from the value zero for the first (few) point (s) of the tra- 
jectory, if the area is sufficiently large the value of df is 
expected to be the value of the stable index, df = a, 
as long as the boundary effects do not come into play. 
In our simulations we do not see a saturation close to 
df = a, as we are interested in the finite size effects and 
therefore choose quite small domains. Due to this finitc- 
ness of the domain the Levy flight will eventually cover 
the entire area and thus the fractal dimension is expected 
to saturate at the value df = 2. This corresponds to the 
observation of ergodicity for bounded Levy flights (2f| • 

To determine the time evolution of the fractal dimen- 
sion of the Levy flight under consideration we calculate 
the box counting dimension of the trajectory [34| . In this 
method, the number of boxes N containing a point of the 
fractal set (here, the points of visitation of our trajectory) 
is counted for a given box size e 2 . For an underlying frac- 
tal geometry and sufficiently small e a power law relation 
of the form N(e) ~ e~ d i is expected [Hj. To validate 
our method we first determined the fractal dimension of 
a number of unbounded Levy flights, for which the value 
df = a is expected. The results are summarized in Tab. HI 
We find that the fractal dimension is reproduced quite re- 
liably. However, the number of necessary random walk 
steps to obtain sufficiently good results increases dramat- 
ically with growing a. 

In Fig. [2j we plotted the time-dependent fractal di- 
mension calculated via the two-dimensional box-counting 
method for confined Levy flights with different values of a 
(a = 0.7, 1.2, 1.7, and 2) for fixed box size 1, 000 x 1, 000. 
As expected for all a values the fractal dimension reaches 
the Euclidean embedding dimension, df — > 2. We see 
that for growing a indeed the slope of the curve df(t) is 
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FIG. 2: Fractal dimension (df) of Levy flights on the square 
domain with periodic boundary conditions as function of time 
for different stable indices a, with a — V2 and a — 500. Here, 
200 sample trajectories were used for the ensemble average to 
produce df at a given time t. The inset magnifies the plateau 
region for a = 2 compared to the case a — 1.7. 

smaller, such that the stationary limit is reached slower. 
Interestingly, for the Gaussian limit a — > 2 an interme- 
diate plateau emerges, as indicated by the arrow in the 
inset of Fig. [2j By analysis of a large number of Levy 
flights for this a value we made sure that this effect is 
not an artifact. The significance of this observation is 
also underlined by the relatively small statistical sample 
error bars shown in Fig. [2] Thus, close to a = 2 care 
should be generally taken when determining the fractal 
graph dimension of a Levy flight process. 

The relaxation time to the ergodic behavior of full area 
coverage of a Levy flight may reach considerable values. 
For a searching animal it is therefore imperative to have 
a reasonable field of vision, i.e., to be able to scan an area 
around each point of the trajectory. Even more efficient is 
constant scanning as discussed in Ref. Q, such that the 
scanned on-the-fly area represents a sausage. Without 
such provisions the scanning of the area by the Levy flight 
remains sparse a forteriori, unless extremely long search 
times arc invested. 



B. Probability density function 

We now turn to the probability density function of 
the process. Fig. [3] shows the time evolution of the 
unidirectional probability density function P(x, t) for 
Levy flights with stable indices a = 0.7 and 1.7. For 
comparison we also present graphs for Brownian mo- 
tion, a = 2. Starting with the sharp initial condition 
Pq(x) = limj-^o P{x, t) = <5(x), the probability density 
function successively broadens, until it reaches the state 
of equidistribution, P s t(x) = l/(2a). We observe that the 
specific characteristics of the different cases < a < 1, 
1 < a < 2, and a = 2 remain discernible up to rela- 
tively long times. Thus up to some 10% of the overall 
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FIG. 3: Time evolution of the unidirectional probability dis- 
tribution function P(x, i) of confined Levy nights with a = 0.7 
and 1.7, as well as Brownian motion with a — 2 in a box of 
edge length, 2a — 1,000, for a — \/2. In each case 10,000 
sample trajectories were used to create an individual curve 
for P(x,t). Fits (dashed black lines) correspond to stable 
distributions as in Eq. lj2]). 



simulations time, it is still possible to distinguish the ba- 
sic shapes of the probability density functions for these 
cases. The overall simulations time needed to reach satu- 
ration is, however, dramatically different, as already ob- 
served in the labels shown in Fig. [3l This Figure also 
demonstrates how fitting to a stable law, Eq. ([2]), yields 
an approximation of a. The fits in this figure show ex- 
cellent agreement with simulation data. 

Thus, from fitting probability density functions alone 
we could conjecture that the initial dynamics of the 




FIG. 4: Tail analysis for the probability distribution function 
P(x, t) of confined Levy flights at an early stage of the random 
walk. The stable fits (dashed black curve) were performed at 
the same time instants as in Fig. [31 but a larger number of 
3 x 10 6 of random walks were needed to produce precise tail 
data (continuous red curve). 



searcher starting far off the edges of the box are essen- 
tially unaffected by boundary conditions: at small times, 
the spatial distribution can be approximated by the sta- 
ble law from an unbounded Levy flight. However, the ef- 
fect of boundaries is known to be highly non-local if the 
random walk is governed by long-tailed jump statistics 
([1]). Due to the non- negligible probability for extremely 
long excursions, the searcher probes the nature of the 
boundaries from the very start of its motion. For exam- 
ple, this leads to modifications on the associated diffusion 
equations dHHH. Fig. @] highlights the specific tail prop- 
erties of the early-time distribution. Despite the almost 
perfect fit quality of stable distributions along virtually 
the whole search space, there are considerable deviations 
in the far tails near the boundaries. In particular, the 
distinct heavy-tail property, Eq. ([T]), is already lost at 
this early stage of the random walk. 

Analyzing the time evolution of moments reveals fur- 
ther distinct peculiarities of the bounded motion, as we 
show in the following. 



C. Moment analysis 

First, consider the mean squared displacement (x 2 (t)), 
which is another common approach to assess the area 
coverage of the process. While this quantity diverges 
for a free Levy flight or a Levy flight in an harmonic 
external potential [H, HH, it is already finite in steeper 
than harmonic external potentials [26| . In our finite area 
the mean squared displacement is thus a valid measure 
for the motion. Figure [5] shows the time dependence of 
the one-dimensional variance, for different values of the 
stable index a. The saturation of the mean squared dis- 
placement is comparatively sharp, such that we can read 
off a typical time scale r from Fig. [SJ sec below. In the 
initial diffusion regime, t <C t, all graphs depict a lin- 
ear time dependence, independently of a. This highly 
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FIG. 5: One-dimensional mean squared displacement of Levy 
flights in a box with periodic boundary conditions and size 
1,000 x 1,000 as function of time t, for different values of 
a. Here, 10,000 sample trajectories have been used for the 
ensemble average {x 2 (t)}. The horizontal black line represents 
the expected stationary value a 2 /3 = 83, 333, where a = 500. 
Note that the early diffusion process is "normal" in the sense 
that (x 2 (t)) grows linearly in time, despite the underlying 
scale-free jump distributions. 
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FIG. 7: Scaling exponent e(S) for moments {\x\ s } of arbitrary 
order 8 > 0. The exponent is obtained by fitting to a power 
law, (\x(t)\ s ) ~t e(l5) . The theoretical prediction (black dotted 
lines, Eq. (O) is valid only in the asymptotic limit 1< t < 
t ~ a°°. For this simulation, 2a = 10 5 , t £ {10, 30}. The 
calculation of higher order moments exploits the far tails of 
the distribution, requiring a relatively large number of 10 8 
trajectories. 
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FIG. 6: Diffusion constant versus box size. With the initially 
"normal" growth of the mean squared displacement, one can 
associate a diffusion constant D x = {x 2 (t)) / (2t) . The latter 
has a power law dependence on the edge length 2a of the 
box. Only for bounded Brownian motion, a — 2, the (large) 
system size does not affect the initial diffusion parameters. 
Error bars are of the size of the symbols. 



contrasts the properties of an unbounded Levy flight, 
where x 2 scales like t 2 l a . In fact, from measuring the 
mean squared displacement alone one could assume the 
motion to be confined Brownian. Thus, the effect of the 
boundaries goes beyond ensuring the finiteness of the sec- 
ond moment. Interestingly, the boundary characteristics 
also have a distinct quantitative effect on the initial dy- 
namics, which we see in Fig. [5] the associated diffusion 



coefficient D x = (x 2 (t))/(2t) depends on the size of the 
system. Enlarging the box increases the (finite) variance 
of individual jump lengths, thereby increasing D x . 

To get the full picture we compute arbitrary mo- 
ments of absolute displacements from our simulations. 
We generally find an initial power law time dependence, 
~ i e W, before the occurrence of saturation. 
More precisely, from the data for the scaling exponent 
e{5) displayed in Fig. we expect that in the asymp- 
totic limit 1 < t < r 

!t 5 / a 7 0<cv<2, < (5 < a, 
t, 0<a<2, 5>a, (3) 
t 5 / 2 , a = 2, S>0. 

We note that this type of moment scaling is actually 
a familiar one in the theory of Levy flights. It also oc- 
curs when truncating the distribution ([1]) for the jump 
lengths at some fixed distance j37H39ll. or when consider- 
ing the effects of a finite ensemble [40||. Roughly speak- 
ing, the truncation distance is analogous to the largest 
value drawn from a finite ensemble; in our system this 
role is played by the finite edge length of the box. In 
all such cases, the initial dynamics resemble an ordinary, 
unbounded Levy flight, albeit with finite moments of the 
type ©. 

In particular, we conclude that the second moment of 
bounded motion mimics "normal" diffusion, but the as- 
sociated diffusion coefficient must depend on the system 
size parameter a: at early times, moments with < 5 < a 
are adopted from the unbounded Levy flight and are thus 
independent of a. At late times, all moments saturate at 
(|x| 5 (oo)) ~ a s . The typical time scale for the interaction 
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with the boundaries therefore scales as 



r~a a , but also r ~ a 2 /(2D x ). (4) 

This necessarily implies that 

D x =D x (a)~a 2 - a . (5) 

The above relation is in nice agreement with our power 
law fits in Fig. [6] The exponent also consistently co- 
incides with the one found for the dependence on the 
truncation parameter in truncated, but unbounded, Levy 
nights as anticipated in Ref. (39| . 

D. Mean first passage time 

We finally analyze the mean first passage time (ti) of 
our confined Levy flights to the boundary of the area. 
On a semi-infinite domain Levy flights show the univer- 
sal Sparre Anderson scaling ~ t x ' of the distribution 
of first passage times, so the mean first passage time 
(ti) diverges [26|, |27|. On our finite domain it converges 
and depends crucially on a. We determined the mean 
first passage time for Levy flights with different values 
of a and for different box sizes (a = 250, 500, 1000 and 
2000). From scaling arguments we would expect that the 
mean first passage time grows like (ti) ~ a a as func- 
tion of the initial distance from the absorbing boundary, 
see Eq. Fig. [5] depicts the rescaled mean first pas- 
sage time a~ a (ri) as function of a. We observe that the 
scaling relation holds provided we can approximate the 
Levy flight by a continuous diffusion process. Conversely, 
by choosing small values for a a we generate a random 
walker that escapes the box after only a few number of 
steps and the universal scaling Q breaks down. From 
the monotonic nature of the graphs we conjecture that - 
for certain fixed parameters a and a - one might find an 
optimal value for the tail parameter a which minimizes 
the mean first passage time. 

We note that the analogous first passage time prob- 
lem in one dimension has been treated elsewhere analyt- 
ically and numerically [4l|, |42j ■ A generalization to our 
two-dimensional problem however is not straightforward: 
Although we fix the scale a = y/2 for individual jump dis- 
tances r , displacements are still affected by random jump 
directions 8. One can argue that in effect, the scale of dis- 
placements along the x—, y— or radial coordinate is not 
fixed, but has a non-trivial relation to the tail parameter 
a. Consequently, such considerations raise the general 
question of how to appropriately fix a jump length scale 
for variable stable exponents in multi-dimensional Levy 
flights. 

IV. CONCLUSIONS 

Area coverage of stochastic processes is an important 
criteria for the efficiency of search or spreading processes. 
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FIG. 8: Rescaled mean first passage time of confined Levy 
flights as function of the stable exponent a, for different box 
sizes. Here, 10 4 -10 5 sample trajectories have been used for 
each data point. 



Here we investigated the time evolution of the time av- 
erage of the area coverage in finite territories of Levy 
flights. As measures we used the fractal dimension, which 
relaxes to two at sufficiently long times, and moments of 
displacements of the Levy flight. Consistently, we found 
that both quantities show a relatively immediate satura- 
tion behavior, such that a clear time scale can be deter- 
mined for reaching full area coverage. However, this time 
scale crucially depends on the value of the stable index 
a: when a approaches the limiting value two, this time 
scale increases significantly. As expected, smaller values 
for a, allowing for longer jump lengths effect more effi- 
cient area coverage and thus search efficiency. In general, 
the number of steps necessary to reach full area coverage 
is relatively large. Thus searchers using this type of strat- 
egy will likely not be able to achieve absolute certainty to 
find a given target randomly. A field of vision, as intro- 
duced in the search literature, will mellow this problem 
and render Levy flights a very efficient search strategy. 

The analysis of the area coverage is completed with the 
study of the equilibration behavior of the probability den- 
sity function of the process, allowing one to distinguish 
different domains of the stable index a (smaller/larger 
than one) up to relatively long times. Moreover the mean 
first passage time T\ scales with the index a and the sys- 
tem size approximately in the scaling fashion (t±) ~ a a . 

From the data analysis point of view, we find that 
studying one single type of above measures is usually not 
enough to identify a given motion as a bounded Levy 
flight. The fractal dimension does have a characteristic 
behavior for long-tailed jump lengths, but quite extensive 
data is necessary to achieve reliable results from a box 
counting algorithm. Although fitting the spatial distri- 
butions by stable laws might agree nicely at early stages 
of the process, their characteristic heavy tails are quickly 
reshaped by interactions with the boundaries. While the 
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mean squared displacement indicates normal diffusion, 
the associated diffusion constant turns out to depend on 
the system size and lower order moments indicate super- 
diffusion. Finally, an analysis based on first passage times 
only yields decisive results if data on a variety of system 
sizes is available. We conclude that usually a combina- 
tion of methods is necessary to unambiguously determine 
the Levy flight nature of a given process. 
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